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ABSTRACT 



The use and measurement of methyl tertiary butyl ether (MTBE) in 
gasolines is increasing in importance to the oil refining industry. 
Remote single strand fiber optic spectroscopy provides unique advantages 
over currently implemented techniques. Such advantages include: 
accurate results in less than 30 seconds, remote distances of over 50 
meters from analyzer to sample point, and no sample preparation. MTBE 
in the range 0-16 volume % may be measured through Second Overtone Near 
Infrared (SONIR) spectrophotometry and multivariate calibrations with an 
associated standard error of prediction of 0.121 volume % MTBE. This 
technique introduces the potential for on-line applications in oil 
refineries and field checking capabilities for regulatory agencies. 



INTRODUCTION 

In recent years the U.S. gasoline pool has been undergoing multiple 
changes. These changes include: higher octane products, lead 
phasedovm, and lower Reid Vapor Pressure (RVP) products. As a result, 
MTBE is gaining widespread appeal with refiners as a gasoline blending 
component for both premium and regular grade gasolines. This is due to 
its high octane and moderate volatility properties without adversely 
impacting drivability. The currently accepted technique for measuring 
MTBE is the "Standard Test Method for Determination of C^ - C 4 Alcohols 
and MTBE in Gasolines by Gas Chromatography" (ASTM D4815T l1 ' , where MTBE 
is determined in the concentration range 0.1 to 10.0 volume percent - 
Precision statistics associated with this method are: reproducibility = 
24% and repeatability = 10% as defined by ASTM criteria 11 ' . NIR fiber 
optic spectrophotometry with multivariate calibration provides a rapid 
and reliable alternative to the gas chromatographic determination of 
MTBE in gasolines. 

Partial least squares (PLS) is the multivariate calibration technique 
used throughout this paper. PLS permits multivariate calibration of NIR 
analyzers by mathematically recording a calibration model. A PLS model 
contains correlation information from two matrices. In this application 
the two matrices are: 1) a set of spectra and 2) a set of corresponding 
MTBE concentrations. Once a model is generated, spectra are used as 
input data in model computations of MTBE. The premise for using SONIR 
data in modeling MTBE, is the spectral bands which are characteristic of 
MTBE. The internal oxygen present in MTBE significantly changes the 
location of the second C-H overtones relative to the C-H overtones from 
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methyl groups attached to primary and tertiary carbons in analog 
compounds that do not contain oxygen. 

EXPERIMENTAL 

Spectra from 74 samples were used in generating PLS models to predict 
volume percentages of MTBE . The spectra were recorded in the range 9 00- 
1600 tun by a double beam scanning fiber optic spectrophotometer 
controlled by a remote host computer (see figure 1} . Two sets of 
spectra (lab and field) were obtained by using 112 meters of single 
strand 500 micron fiber optic cable with a total of three coupling 
points post instrument. One point of fiber coupling occurred at the 
instrument/ sample interface which was accomplished through a 
transmission probe assembly specifically designed to utilize 1 cm 
cuvettes (see figure 2) . The other two points of coupling occurred in 
fiber runs: one in the transmitting and one in the receiving fiber. 
Spectra were recorded against an air reference using 4 scan averaging at 
23°C. A total of 226 wavelengths (1100-1550 nm in even increments) were 
selected for PLS modeling of MTBE. The 1100-1300 nm region contains the 
second overtones from symmetric and asymmetric methyl, methylene, 
methine, aromatic, and olefinic C-H stretches. The 1350-1500 nm region 
contains absorption bands which result from combinations of stretching 
and bending modes from methylene, methine, aromatic, and olefinic C-H 
bonds . Data pre-treatment consisted of two 17 point Savitzky-Golay 
smooths and baseline correction at 970 nm. This method is very similar 
to that used in the NIR determination of octane number in gasolines by 
PLS 13 ' . 

LABORATORY SAMPLES FOR PLS CALIBRATION 

1, 2, 4 , 9, 12, and 16 volume % MTBE additions were made to three 
grades (84, 87, and 91 R+M/2) of unleaded gasoline from a single vendor. 
These 21 (18 with MTBE added and 3 neat) laboratory samples were used in 
the preliminary evaluation of the relationship between the SONIR spectra 
and the volume % MTBE present in each fuel. The SEE from a 4 factor PLS 
model is 0.092 volume % MTBE calculated from (1) . 

SEE = (y jfa - y j/p ) 2 /(n-l)l 1/2 (1) 

j=l 

where: n = number of samples in calibration, 
j * sample number 
yj a = actual value. 
7j p- predicted value. 



An estimate of the SEP (see equation 2) from this modeling was obtained 
by validation. This was accomplished by removing a minimum of one 
sample and a maximum of three from the data set, and generating eight 4 
factor models from the balance of spectra in each case. From this set 
of 21 prediction results, the SEP is 0.121 volume % MTBE. 



SEP * <y j/a - y^p) 2 /^-!)] 1 ' 2 (2) 

j = l 

where: n = number of samples in predictions, 
j = sample number 
y.i a = actual value, 
yj' = predicted value. 
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Figure 3 gives a plot of the actual volume % MTBE against the predicted 
volume % MTBE used to determine the SEP from the laboratory samples. 
This laboratory work provides an estimate of the errors associated with 
instrumental/PLS system. The ability to effectively determine MTBE in 
gasolines using laboratory samples motivated further investigation of 
PI^ modeling using spectra from field samples and volume % MTBE obtained 
by current ASTM methodology* 1 - 



FIEXJ3 SAMPLES 




the ASTM method, where the volume % MTBE was reported as 10.1 and 11.2. 
These results are within and typical of the ASTM repeatability 
statistics from the gas chromatographic method. Both samples and their 
corresponding reported values were included in subsequent calibrations. 
The SONIR prediction results for the duplicate samples are 10.0 and 
10 7 The trends from the gas chromatographic and SONIR analysis using 
the duplicate samples are similar. Since the SONIR results from the 
duplicate samples produced a deviation far greater than expected from 
th% SEP determined using the lab samples, there is an indication that the 
two samples may not be true duplicates. 

The 21 samples from the laboratory set were added to the field samples 
to qive a total of 74 samples for evaluation. The SEE from a 6 factor 
PLS model is 0.8 volume % MTBE calculated from (1). An estimate of the 
SEP (see equation 2) from this modeling was obtained by removing a 
minimum of 6 samples and a maximum of 14 from the data set, an °W om <_ hfa 
generating eight 6 factor models from the balance of spectra. From the 
let o£ 74 prediction results, the SEP is 1.1 volume % MTBE. Figure 4 
gives a plot of the actual volume % MTBE against the predicted volume % 
MTBE used to determine the SEP from the field samples. Due to the 
distribution of MTBE concentrations normally encountered m the fieia 
(either in the 12 % range or no MTBE), a sub-set was selected for 
calibration where the set of samples containing no MTBE (called the nil 
letPwas removed. An estimate of the SEP from the sub-set ""djljn* w « & 
obtained by removing a minimum of 7 samples and a maximum of 14 from the 
data set, and generating five 6 factor models from the balance of 
spectra? From the set of 56 prediction results, the SEP is 0.7 volume % 
MTBE. Figure 5 gives a plot of the actual volume % MTBE against the 
predicted volume % MTBE from predictions using this sub-set. This was 
done in an attempt to illustrate the modeling «P* b i lxt i e ? °f * h *| data 
method without being limited by a poorly distributed population of data. 

A summary of the results obtained from the three sets evaluated is shown 
in table 1. 



Table 1 

Comparison of Results Using Various Data Sets 



Data Set Factors Used SEE SEP . 

0.092 0.121 



Lab 4 

Field + Lab # 6 0.8 1-1 

Field + Lab* 6 0.6 0.7 

* The nil set from the field + lab samples were removed for this 
analysis. 

The total number of factors used to predict volume % MTBE increased from 
4 to 6 when field samples were included. This was due £ fc he xncreaae 
in compositional variations of the base gasolines from the different 
vendors (as indicated by the scores statistics from the PLS model) . 
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Notice that the errors in the predictions from the field set are larger 
than the errors in the predictions from the laboratory set. This is 
because the errors in the data from the field reference method (by gas 
chromatography} exceed those of the data from the laboratory reference 
method {by gravimetric addition) - The errors associated with the 
concentrations of the standards used for calibration ultimately limit 
the prediction results from modeling of such referee data. 



THE EFFECTS OF VARIABLE FIBER LENGTHS 

The fiber optics used in this experiment exhibited a silanol absorption 
band* 5 ' at 138 5 nm of approximately 1 absorbance unit per 100-meters of 
fiber. One factor which should be addressed when using fiber optics is 
the effect of such absorption bands on the measured spectrum. In this 
application of fiber optic spectrophotometry, the silanol band occurs in 
the range of wavelengths used for calibration. One effect of this 
absorption band is the addition of photometric noise due to the decrease 
in absolute light levels reaching the detector. Since the 112 meters of 
fiber previously used gave a silanol band with a peak absorbance of 
approximately 1.1 absorbance units, a shorter run of fiber was selected 
in an attempt to reduce the peak absorbance of the silanol band. This 
should improve the signal-to-noise characteristics of the spectra 
exclusively in the silanol band region, and possibly improve the 
modeling of MTBE. The samples from the laboratory calibrations were 
used to determine the SEE and SEP from calibrations with spectra 
recorded by using 12 meters of fiber optic cable. 

Figure 6 gives a plot of the actual volume % MTBE against the predicted 
volume % MTBE from the work done using 12 meters of fiber. 

Table 2 gives a summary of the comparison results obtained using 12 and 
112 meters of fiber optic cable. 



Table 2 

Comparison of Results from Different Fiber Distances 

Fiber Distance, m SEE SEP 

12 0 .092 0.126 

112 0.108 0.121 

The 112 meters of fiber optic cable used in this comparison does not 
appear to limit the modeling of MTBE. 

CONCLUSIONS 

The statistical results obtained from independently modeling volumetric 
and gas chromatographic data from several sets of samples indicates that 
the limiting factor is the uncertainty associated with the currently 
implemented gas chromatographic method- SONIR spectroscopy with 
mutivariate calibration techniques introduces an accurate and rapid 
method of monitoring levels of MTBE for production and laboratory 
purposes. This coupled with the use of long fiber lengths provides 
remote, intrinsically safe installation of the NIR analyzer (such as 
control room location) , which permits convenient on-line monitoring of 
MTBE for refineries. 

An additional benifit of remote SONIR spectroscopy which is often 
overlooked is the elimination of sampling problems- In this study a 
potential sampling problem was highlighted by the gas chromatographic 
and SONIR results obtained using the duplicate samples- The samples 
used to evaluate the repeatability of the ASTM and SONIR methods appear 
to be individuals rather than duplicates. Since sampling is not 
required with an operational on-line measurement system, this problem 
would be eliminated. Furthermore, this technique also has potential for 
regulatory agencies , where the equipment similar to that used in this 
experiment could be modified for mobile, field monitoring applications - 
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